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Abstract. 

We provide conditions for a lattice scheme defined on a four points lattice to be 
linearizable by a point transformation. We apply the obtained conditions to a sym- 
metry preserving difference scheme for the potential Burgers introduced by Dorod- 
nitsyn 1 3 1 and show that it is not linearizable. 
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1. Introduction 

In a recent article |@] we extended to lattice equations the theorems introduced by 
Bluman and Kumei [2] for proving the linearizability of nonlinear Partial Differen- 
tial Equations (PDEs) (for a recent extended review see 03) based on the analysis 
of the symmetry properties of linear PDEs. In [£[], following the analogy of the 
continuous case we formulated a similar theory for linearizable Partial Difference 
Equations (PAEs) on a fixed lattice. 

Here we extend the results of [4 1 to the case of a lattice scheme, i.e. when the lattice 
is not a priory given but it is defined by an equation so as to be able to perform a 
symmetry preserving discretization of the PDE. 

In Section 2 we prove a theorem characterizing the symmetries of linear difference 
schemes on four lattice points and in Section 3 we apply it to find conditions under 
which a nonlinear difference scheme is linearizable. These conditions are then 
applied to the symmetry preserving discretization of the potential Burgers. 

2. Symmetries of linear schemes 

In this Section we define a difference scheme and provide the symmetry conditions 
under which such a scheme is linearizable. To do so in a definite way we limit 
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ourselves to the case when the equation and the lattice are defined on four lattice 
points in the plane, i.e. we consider one scalar equation for a continuous function 
of two (continuous) variables: u m n = n(x m ra , t m n ) defined on four lattice points. 



Um+l,n,m 



Figure 1. The 1? square-lattice where the equation is defined 



2.1. The difference scheme 



As we consider one scalar equation for a continuous function of two (continuous) 
variables, a lattice will be a set of points Pi, lying in the plane M? and stretching 
in all directions with, a priori, no boundaries. The points Pi in M 2 will be labeled 
by two discrete labels P m ,n- The Cartesian coordinates of the point P m ,n will be 
(x m ,n,t m ,n) with — oo < m < oo , — oo < n < oo. The value of the dependent 
variable in the point P m>n will be denoted u m>n = u{x m ^ n , i m , n )- 
A difference scheme will be a set of b equations relating the values of {x, t, u} in 
a finite number of points. We start with one 'reference point' P m ,n and define a 
finite number of points 

Pm+i,n+j hi the neighborhood of P m ,n- 

They must lie on 

two different sets of curves, two of which will be intersecting in P m ,n- Thus, the 
difference scheme will have the form 



{ x m+i,n+ji tm+i,n+ji u m+i,n+j} ^ — 1 < a < 6 

-h<i<i 2 -ji<j<32 n,«2,Ji,j2 e 



(2.1) 



The situation is illustrated on Figure 2 in the case of a present 7 points lattice. Our 
convention is that x increases as m grows, t increases as n grows (i.e. x m+ i :Tl — 
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x m ,n = hi > , t mn+ i — t mn = h 2 > 0). The scheme on Figure 2 could be 
used e.g. to approximate a differential equation of third order in x, second in t. 



The value of b, the maximum number of different equations we consider, depends 
on the kind of problems we are considering. Starting from the reference point 
P m ,n and a given number of neighboring points, it must be possible to calculate 
the values of {x, t, u] in all points in a unique way. This requires a minimum 
of three equations to calculate the independent variables (x, t) in two directions 
and the dependent variable u in all points. With one dependent variable in M. 2 , at 
most we can set b = 5. Of the five equations in (12- lb . four determine completely 
the lattice, one the difference equation. If we choose b = 3 than two define the 
lattice and one the difference equation and we are solving an initial value problem 
when both the equation and the lattice are defined from given initial conditions. If 
a continuous limit exists, (12.11) represent a PDE in two variables. The equations 
determining the lattice will reduce to identities (like = 0). 

As an example of difference scheme, let us consider the simplest and most standard 
lattice, namely a uniformly spaced orthogonal lattice and a difference equation 
approximating the linear heat equation on this lattice. The 5 equations (12.11 ) in this 
case are: 



t 




m+2,n 
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Figure 2. Points on a two dimensional lattice 
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Figure 3. Four points in the case of the heat equation. 



(2.3) 



Um,n+1 Umn "^m+l.n 2u mn + Um— l,n 
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(2.4) 



where /ii and /12 are two constants. 

The example is simple; the lattice equations can be solved explicitly to give 



him + xq t r 



h 2 n + 1 . 



(2.5) 



The usual choice is xq = to = , hi = h% = 1 and then x is simply identified 
with m, t with n. The example suffices to bring out several points: 

1 . Four equations are needed to describe completely the lattice but in this case 
there is a compatibility condition. In the whole generality two equations are 
sufficient and provide the lattice starting from some initial conditions. 

2. Four points are needed for equations of second order in x, first in t. Only 
three figure in the lattice equation, namely P m +i,n>Pm,n an d Pm,n+i- To 
get the fourth point, P TO _x, n , we shift m down by one unit the equations 
(I2.2H2.41) . 

3. An independence condition is needed to be able to solve for x m+ i n , t m +i, n , 

^m,7i+l> tm,n+l and U m rl -)_i. 



We need the more complicated two index notation to describe arbitrary lattices and 
to formulate the symmetry algorithm. 
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2.2. Symmetries of the difference scheme 

We are interested in point transformations of the type 

x = F\(x,t,u) t = G\(x,t,u) u = H\(x,t,u) (2.6) 



where A is a group parameter, such that when (x,t,u) satisfy the system (12.11 ) 
then (x, t, u) satisfy the same system. The transformation acts on the entire space 
(x,t, u), at least locally, i.e. in some neighborhood of the reference point P m ,n, 
including all points P m+ j „ + j figuring in equation (I2.1I ). That means that the same 
functions F, G and H determine the transformation of all points. The transforma- 
tions ( 12.61 ) are generated by the vector field 

X = £(x, t, u)d x + t(x, t, u)d t + <j>(x, t, u)d u . (2.7) 

The symmetry algebra of the system (12.11 ) is the Lie algebra of the local symmetry 
group of local point transformations. An infinitesimal symmetry (12.71 ) is a symme- 
try of (12.11 ) if (12.11 ) is invariant under a transformation (12- 6b - To check it we must 
prolong the action of the vector field X from the reference point {x m ^ n ,t m ^ n , u m>n ) 
to all points figuring in the system (12.11 . Since the transformations are given by 
the same functions F, G and H at all points, the prolongation of the vector field 
(12.71 ) is obtained simply by evaluating the functions £, r and cj> at the corresponding 
points. In other words, we have 

(2.8) 



where the summation is over all points figuring in the system (12.11) . The invariance 
requirement is formulated in terms of the prolonged vector field as 

prXE a \ Ec=Q = l<a,c<6. (2.9) 

Just as in the case of PDE's [6], we can turn equation ( 12.91 ) into an algorithm for 
determining the symmetries, i.e. finding the coefficients in vector field (12.71 ) [5]. 



2.3. Symmetries of a linear partial difference scheme 

To be able to linearize a difference scheme using the knowledge of its symmetries 
we must be able to characterize the symmetries of a linear scheme. To do so in 
this subsection we prove a theorem on the structure of the symmetries of a linear 
partial difference scheme: 
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Theorem 2.1 Necessary and sufficient conditions for three difference equations 
&m,n = 0, 3vn, n = and & m ,n = defined on four points {(m, n),(m + 
1, n), (m, n + 1), (to + 1, n + l)}/or a scalar function u m) „(x m) „, * m .„) a«<i the 
lattice variables x m , n and t m<n to be linear is that they are invariant with respect 
to the following infinitesimal generator 

X. m ,n = Vm tn d Um n + Xm,nd Xm ^ n + rjm,ndt m n , (2-10) 

where the discrete functions v m ,n, Xm,n, f]m,n satisfy three linear equations 
v m +l,n+l = e m ,n, Xm+l,n+l = fm,n and T m+ i )n+ i = g m , n . The functions t, f and 
q depend just on the functions (v m ^ n , Xm,n> Vm,n) m the points (m, n), (to + 1, n) 
and (to, n + 1) and are given by: 

eo,o = ait;o,o + 02^0,1 + ^1,0 + 04X0,0 + 05X0,1 + ^6X1,0 + a7%,o+ 

+ 08^0,1 + agr/1,0, 
fo,o = Mo,o + Mo,i + 63^1,0 + 04X0,0 + 05X0,1 + 06X1,0 + Mo,o+ 

+ o 8 %,l + o 9 ??i,o, 

00,0 = civ 0fi + c 2 v 0: i + 03^1,0 + c 4 xo,o + c 5 Xo,i + c 6 Xi,o + c 7 r?o,o+ (2- 1 1) 
+ csVo,i + cg??i,o, 

where a\, ■ ■ ■ , eg depend only on the lattice indices and where, here and in the 
following, for the sake of simplicity we set in any discrete variable on the square 
z m +i, n +j = Zi,j- The linear equations (£ m , n = 0, 3m, n = and & m ,n = have 
the form: 

""1,1 = ai^o.o + a2-"o,i + 03^1,0 + a 4 a;o,o + a 5^o,i + a6^i,o + «7*o,o+ 

+ (28*0,1 + 0-9*1,0) 

= 6i-u ,o + o 2 iio,i + b 3 u lfi + 04x0,0 + 05x0,1 + 6 6 xi, + 67*0,0+ 

+ 08*0,1 + 6g*i,o, 

*i,i = ci-u ,o + c 2 n ,i + c 3 ni, + c 4 x ,o + c 5 x ,i + c 6 xi, + c 7 t ,o+ (2.12) 
+ C8*o,l + Cg*i,o- 

Proof. To prove this Theorem we require that a generic PAE F m n = 0, depend- 
ing on a scalar function u m , n (x m , n , i m , n ) and the lattice variables x m , n and t m n 
in the four points {(m, n), (m + 1, n), (m, n + 1), (m + 1, n + 1)}, i.e. 12 vari- 
ables, be invariant under the prolongation of (12.101 . as given by (I2.8I ). The in- 
variance condition (|Z9]>, when £ m>n (x,t,u) = Xm,n, T m>n (x,t,u) = ?? m , n and 
4'rn,n{x,t,u) = u m n implies that F m n should depend on a set of 11 invariants 
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v m ,n, Xm,n and r? m>n independent 

Li = v ,ou ,i - v ,iu 0! o, L 2 = v fiu li0 - t>i,o«o,o, 

L3 = v ,oUi,i ~ eo,o^o,o, L A = v 0fi x 0jl - Xo,i«o,o, 

L5 = voflxifl - xi,o«o,o, L& = vojox\ t \ - fo,o«o,o, 

L7 = «o,o*o,i - %,i«o,o, ^8 = ^0,0*1,0 - "i,o«o,o, 

^9 = ^0,0*1,1 - 00,0^0,0, -^10 = ^0,0^0,0 — Xo,o«o,o, 

^11 = ^0,0*0,0 - Vo,ouo,o- (2-13) 

As F m „ should not depend on the functions (v mjn , Xm,n, Vm,n) i n the points 
(m, n), (m+1, n) and (m, n+1) we have nine constraints given by » l9 - Fm - n — = 0, 

and .,„fW* . = with = (0, 0), (0, 1), (1, 0). These are 



first order partial differential equations for the function F m n with respect to the 1 1 
invariants which we can solve on the characteristics to define three invariants 

K\ = VQfi{ui : i — [eo,o,u ,i u o,i + eo,o,Di,o u i,o + ^o,o,xo,o x o,o + «o,o,xo,i x o,i+ 
+ e o,o,xi, 0^1,0 + ^0,0,770,0*0,0 + eo,o,»?o,i*o,i + eo,o,??i,o*i,o]} - uo,oUo,o— 

- [60,0,1)0,1^0,1 + eo,o,Di,o^i,o + eo,o,xo,oXo,o + eo,o,xo,iXo,i + eo,o,xi, Xi,o+ 

+ e O,0,7)o, 0^0,0 + e O,0,7)o, 1^0,1 + e 0, 0,771,0 
K 2 = V 0fi {ui : i - [fo,0,« ,i n O,l + fo,0,Di, ^l,0 + fo,0,xo,o x O,0 + fo,0,xo,i x O,l + 

+ fo,o,xi,o x i,o + f 0,0,^,0*0,0 + fo,o,??o,i*o,i + fo,o,m,o*i,o]} - «o,o{fo,o- 

- ^0,0,1-0,1^0,1 + fo,0,t)i,o^l,0 + fo,0,xo,oXO,0 + fo,0,xo,iXO,l + fo,0,xi,oXl,0+ 
+ foAw.o^.o + foAJJo.i^.l + fo,0,7ji,o»7l,o]}, 

#3 = f0,o{«l,l - [0O,O,«o,i n 0,l + 0O,O,di,o«1,O + 0O,O,xo,o x O,O + 0O,O,xo,i a: O,l + 

+ 0o,o,xi,o x i,o + 0o,o,r)o,o*o,o + 0o,o,r?o,i*o,i + 0o,o,m,o*i,o]} - ""o,o{0o,o- 

- [00,0,1-0,1^0,1 + 0O,O,t)i,o u l,O + 0O,O,xo,oXO,O + 0O,O,xo, 1X0,1 + 0O,O,xi, 0X1,0 + 
+ 0O,O,J7o,o ^0,0 + 00,0,770,1 V0,1 + 00,0,7)1,0^1,0]}- (2.14) 

By construction the three invariants Ki, i = 1, 2, 3 are independent and the three 
equations <£ m , n = 0, SVn.n = and m , n = must be defined in terms of them. 
The three invariants K3, L3 and M3 still depend on the functions (v mjJl , Xm,n, 
%i,n) in the points (m, n), (m + 1, n) and (m, n + 1) while they should depend 
just on the variables (u m ^ n , x m . n , i m , n ) in the points (m, n), (m+1, n), (m,n + l) 
and (m + 1, n + 1). The derivatives F mi „ z = 1, 2, 3 will satisfy a set of nine 
linear equations whose coefficients will form a matrix 21 9x3. The matrix 21 can 
have rank 3, 2 or 1. In the case of rank 3 we have F m ^ n ^Ki = 0, i = 1, 2, 3 i.e. the 
function F m n does not depend on the 3 invariants. If the rank of 21 is 2 or 1 we can 
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have at most two independent invariants. If we want to have three invariants we 
need to require that the coefficients of the matrix 21 be zero, i.e. defining a\ = t>o,o> 

dK 

a 2 = v ,i, a 3 = v lfi , • • • , a 9 = r/ lj0 we have = Op = 1, 2, 3, q = 1, • • • , 9. 

The equations -q^- = are linear homogeneous expressions in Uij, Xjj and ij j 
with coefficients depending on Vi j, Xi,j an d for appropriate values of i and j. 
Consequently ( I2.12I ). Than = turn out to be a set of 159 overdetermined 
partial differential equations for the functions e m>n , f m>n an d Qm,n whose solution 
( 12.111 ) is obtained using Maple. It depends on 27 integration constants which must 
be set equal zero if (12.121 ) does not depend on Vij, Xi,j an d Tfrj. I 
A few remarks can be derived from Theorem 12. H and must be stressed. 

Remark 1 The equation for u m ^ n and those for the lattice variables x m ^ n and t m ^ n 
are independent, however the functions appearing in the symmetry (12.101) do not 
satisfy equations independent from those satisfied by the lattice scheme. In fact 
these symmetries correspond to independent superposition laws for the equation 
and the lattice. 

Remark 2 If the linear equation for u mn is autonomous than the coefficients 
{04, ■ ■ ■ , ag} are zero. The variable w m>n will satisfy a similar equation but the 
lattice equations can depend linearly on u m ^ n . 

Remark 3 The proof of Theorem 12. 1 1 does not depends on the position of the four 
lattice points considered, i.e. {(m, n), (m + 1, n), (m, n + 1), {m + 1, n + 1)}. 
The same result is also valid if the four points are put on the triangle shown in Fig. 
3, i.e. {(m, n), (m + 1, n), (m — 1, n), (to, n + 1)}. 

3. Linearizable nonlinear schemes 

In this article each equation of a difference scheme is an equation for the contin- 
uous variable u m ^ n , x m ^ n and t mj „. If the equations for the lattice variables, x m n 
and t m)n , are solvable we get 

x m , n = X(m, n,co,ci, ■■■ ), t m>n = T(m,n,d 0) di,---)> (3.1) 

and then the remaining equation for the variable u m ^ n depends explicitly on n and 
to, on the integration constants contained in (13.11 ) and turns out to be an algebraic, 
maybe transcendental, equation of u m . n in the various lattice points involved in 
the equation. So the difference scheme reduce to a non autonomous equation on a 
fixed lattice and for its linearization we can apply the results of [4J. 
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If the equations for the lattice are not solvable the difference scheme can be thought 
as a system of coupled equations for the variables u m>n , x m>n and t m>n on a fixed 
lattice. In this way we can apply to the equations of the scheme the results of 0] 
and, taking into account the results of the previous section, we can propose the 
following linearizability theorem: 

Theorem 3.1 A nonlinear difference scheme d2.il ) involving i\ + i2 different points 
in the m index and j\ + j'2 in the n index for a scalar function u m ^ n of a 2— 
dimensional space of coordinates x m n and t m , n will be linearizable by a point 
transformation 

W m ,niym,ni ^m.n) — f { 1 Xm,n^m,n^m,n) i Vm,n — 9{Xm,ni ^m,n> ^m,n) > (3-2) 

to a linear difference scheme of the kind of \2.12\ for w m ^ n , y m ,n an d z m ^ n if it 
possesses a symmetry generator 

X = £(x,t,u)d x + 4>(x,t,u)d t +ifj(x,t,u)d u , (3.3) 
£(x,u) = a(x,t,u)y, cp(x,t,u) = (3(x,t,u)z, tp(x,t,u) = j(x,t,u)w 

with a, (3 and 7 given functions of their arguments and y, z and w an arbitrary 
solution of <\2.m . 

In the following we will consider the application of this theorem to a difference 
scheme which one would hope that it is linearizable as is a symmetry preserving 
discretization of a linearizable PDE, the potential Burgers equation O. 



3.1. Application 



We consider here the discretization of the potential Burgers presented by Dorodnit- 
syn et. al. |3 ] and show that, even if it is reducible by a point transformation to the 
discrete scheme of the heat equation, is not linearizable by a point transformation. 
As a consequence we have that also the symmetry preserving discretization of the 
heat equation presented by Dorodnitsyn et. al. is not a linear difference scheme. 

The symmetry preserving discretization of the potential Burgers is given by the 
following scheme 



Ax 

T 



t"m,n+l 



1 



h+ + h- 



7? 



w) + 



(W — W- 



w + 



w 



w — w. 



{h+f L 



h + 
t"in—l,n 



(3.4) 

(3.5) 
(3.6) 
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where r is a constant and 



W 



W 



m+l,n j 



Ax — X m ,n+1 %m,ni h — X m +l,n %m,ni h — %m,n %m—l,n- 

Eqs. ( 13.41 13.51 ) are written in terms of the discrete invariants X 2 , Z 3 , Z 4 on the 
stencil defined in terms of (r, x, Ax, h + , h~, w, w, u>+, W-) of the finite point 
symmetries of continuous potential Burgers equation 

1 



wt = w xx 



2^' 



Xi = d t , X 2 = d x , X 3 = td x + xd t , X 4 = 2td t + xd x , 



(3.7) 
(3.8) 



X5 = d u 



1-1 — "TIT; ^2 
h 

Is = ~. + 

4 r 

Z4 = Ax — 
r 



X 6 = t 2 d t + txd x + (^x 2 + / ) i),,.. 

= I e |(t«-«;)+ V 

/l+ 2 



e2 l 



(3.9) 



2h+ 



w 



W — W- 



h+ h- 

— (w + -w) + —{w - W-) 



h+ + h- 

and goes into it in the continuos limit. 

Eqs. (I3.4[ 13.51) are related to a symmetry preserving discretization of the heat 
equation for u m ^ n by the point transformation w mjTl = — 21og(it m>n ). However it 
is not completely obvious if d3.4M3.5b are reducible to a linear discrete equation i.e. 
if it possess, as its continuos counterpart (I3.71 ). an infinite dimensional symmetry 
X = u(x, t)e~ w d w with u(x, t) solution of the linear heat equation ut = u xx . 

We can apply on the lattice scheme (13.41 13-51 13.61 ) the symmetry generator 

X = ip(x, t, w)ud w + 4>(x, t, w)sdt + £(x, t, w)yd x , (3.10) 

with (x,t,w) satisfying (13.41 13.51 13.61 ) while (y,s,u) are solutions of the linear 
scheme prescribed by Theorem 12. II 

U m ,n+1 = aiU m ,n + 0-2U m -\,n + a3«m+l,n + 0>Aym,n + 0-hVm-\,n + fl6?/m+l,ra 
+0-7 s m,n + a8S m _l jn + agS m +i jn , (3.11) 
ym,n+l = ClUm,n Q^m— l,n + C^U m -^i^ n + C/±y m ^ n + C^y m —\^ n 

+ Mm-l,n + ^3 u m+l,n + &42/m,n + &52/to-1,tj + ^62/m+l,n 
+&7S °m-l,n 
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where (a,-, bj, cj, j = 1, • • • ,9) are parameters at most depending on n and m. 
By a long and tedious calculation carried out using a symbolic calculation program 
we get that 

tl>(x,t,w) = Mt) + Mt)x + Mt)x 2 , (3.12) 

fax,t,w) = 4> (t) + 4>i(t)x + 4> 2 (t)x 2 , 
£(x,t,w)=£ (t)+£ 1 (t)x. 

Introducing (13.121 ) into the determining equations for the symmetries of the discrete 
potential Burgers scheme (13.41 13.51 13.61 ) we get 1672 equations for the functions 

(ijjj(t) : 4>j{t), £j(t), j = 0,1,2) depending on the coefficients (aj, bj, Cj, j = 
1, • • • ,9). 168 of those equations do not depend on the coefficients (a,j, bj, cu, 
j = 1, ■ ■ ■ ,9) and on (ipj{t + t), cf)j(t + r), £j(t + r), j = 0, 1, 2); solving them 
imposing that r / Owe get ipj(t) = for j = 0, 1, 2, fa = for k = 1, 2 and 
£fc = for k = 0, 1. Introducing this result in the remaining 1508 equations, we 
get the following 9 equations 

bi<t> (t + T ) = hMt + r) = 6 3 0o(i + r) = & 4 <fo(* + T ) = hM* + t) = 
= beM* + t) = 0o (t) - &70o(i + r) = 6 8 0o(i + r) = 6 9 0o(i + r) = 0. 

If we require (j)o(t) be not identically null, the coefficients bj, j = 1, • • • 6, 8, 9 
must be all zero and b-j 7^ 0. As a consequence 0o(t) = b^ n (j), with (f> an arbitrary 
constant. In this case we have a symmetry generator X = b^ m sdt which is a 
consequence of the linearity of (I3.6I ). So we can conclude that the potential Burgers 
scheme (13.41 13.51) is not linearizable and that the corresponding discretization of 
the heat equation Q is not given by a linear scheme. The linearity of the lattice 
equation for t mn (13.61) is confirmed by the presence of the symmetry X = bj n sdt. 
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